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THE IMPACT OF THE CONVECTIVE BLUESHIFT EFFECT 
ON SPECTROSCOPIC PLANETARY TRANSITS 

Avi Shporer^'^, Tim Brown^'^ 
ABSTRACT 

We present here a small anomalous radial velocity (RV) signal expected to be 
present in RV curves measured during planetary transits. This signal is induced 
by the convective blueshift (CB) effect — a net blueshift emanating from the 
stellar surface, resulting from a larger contribution of rising hot and bright gas 
relative to the colder and darker sinking gas. Since the CB radial component 
varies across the stellar surface, the light blocked by the planet during a transit 
will have a varying RV component, resulting in a small shift of the measured 
RVs. The CB-induced anomalous RV curve is different than, and independent 
of, the well known Rossiter-McLaughlin (RM) effect, where the latter is used for 
determining the sky-projected angle between the host star rotation axis and the 
planet's orbital angular momentum axis. The observed RV curve is the sum of 
the CB and RM signals, and they are both superposed on the orbital Keplerian 
curve. If not accounted for, the presence of the CB RV signal in the spectroscopic 
transit RV curve may bias the estimate of the spin-orbit angle. In addition, future 
very high precision RVs will allow the use of transiting planets to study the CB 
of their host stars. 

Subject headings: planetary systems 



1. Introduction 

Radial velocity (RV) measurements of Sun-like stars are now reaching th e 1 m s~^ accu 



racy milestone and are expected to become more accurate in the future (e.g.. iLi et al.ll2008 
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Dumusque et al.l I2OIOI : iHoward et al.l I2OIOI : ILo Curto et al.l |2010[ ) . Naturally, the increased 
sensitivity allows to identify more and more low amplitude effects that need to be accounted 
for when modeling the measured RV curves. We present here such an effect, induced by 
what is known as stellar convective blueshift (hereafter CB; explained in more detail below), 
which manifests itself as a small anomalous RV signal during a planetary transit. 

This small RV signal is super posed on, and independent of, th e already well known 



Rossiter-McLaughlin (RM) effect (IRossiter 



seen during a planetary transit (e.g., iQueloz et al. 



1924 



2000 



McLaughhnI Il924l: see also iHoltl ll893|) 



Winn et al.ll2005f l. The RM effect 



originates from the host star rotation, and allows the measurement of A, the sky-projected 
angle between the host star's rotation axis and the planet's orbital angular momentum^, 
which is an im portant piece of inforniation i n the study of planetary formation a nd orbital 



evolution (e.g.. iFabrycky fc Tremaind 120071 : iTriaud et al.ll2010l : IWinn et al.ll2010l ). During 



transit the planet blocks the light coming from different regions of the host star surface, with 
a different RV component due to the star's rotation. The observed line profile is distorted, it 
becomes asymmetric, and its center is shifted, resulting in the RM RV signal. The RM RV 
curve shape depends primarily on the spin-orbit a ngle A, the sky-project ed rotation velocity 
Vroti and the planet to star radii ratio r = Rp/Rs (iGaudi fc WinnI 120071 ). 

The new effect we present here is expected to be at the 1 m s~^ level and to occur 
during transit, thus affecting the interpretation of the RV curve. We describe the CB effect 
in Section m and in Section Owe use a simple numerical model to calculate the impact of the 
CB effect on RVs measured during transit. We give a discussion in Section [Hand a summary 
in Section [51 



2. The Convective Blueshift (CB) Effect 



The history of the convective blueshift effect dates back to when astronomers started 
investigating in detail the spectrum of the Sun, and identified disagreements betwe en the 
meas ured hne positio ns and their expected positions based on lab experiments (jjewell 
18961 : for a review see lDravinslll982l ). We know today that several physical processes con- 
tribute to these disagreern ents, one of them being the relativistic gravitational redshift (e.g., 
Earman fc Glymourlll980l ). Still, even after the gravitational redshift and the Doppler shift 
induced by the rotation of the Sun and the Earth were accounted for the above discrepancy 
was not completely resolved. The corrected spectra showed that spectra of the Sun observed 
close to the Sun's limb are redshifted relative to spectra from the center of the Sun's disk. 



Some authors use (3 = 



—A to mark this angle. 
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and this redshift scales with fi = cos 6, where 6 is the angle between the normal to the Sun's 



surface and the line of sight. This u nexpected trend was studied by many (e.g., lEvershed 
19361 : lForbeslll962l : I Adam et al.lll976l ) and it s origin was a puzzle for the better part of the 
20th century (for an early historic review see lForbeslll96ll ). Before it was understood it was 
named the solar redshift problem, or the solar limb effect. 

The solution to the problem came from the realization that the surface of a Sun-like star 
is composed of convective cells, or granules, where hot gas is rising upward, and cooler gas 
is sinking downward in the intragranular lanes. Therefore, the stellar surface is composed 
of local blueshifts and redshifts. Since the rising stream of hot gas is br i ghter than th e 



cooler sinking gas, a net blueshift is observed (e.g., iBeckers fc Nelson! 1 19781 : lDravinslll982l ). 
Following the understanding of the physical process responsible for the effect it was referred 
to as convective blueshift. We denote the "local" blueshift velocity as Vqb, which is the 
blueshift observed at the center of the star's disk, when averaging across a surface area larger 
than the size of a granulation cell. We note that the Vcb velocity is negative. Depending on 
the value of this velocity, the spectral lines are broadene d in an asymm etric way and their 



center is shifted. For the Sun Vcb,© ~ —300 m s ^ (e.g.. lDravinslll987). and it ranges from 



about —200 m s ^ for K-type star s to —1000 m s ^ for F-type stars (jPravins fc Nordlund 



1990al : iDravins fc Nordlundlll990bl : iDravine Il999f). For comparison th e Solar gravitational 
redshift velocity at infinity is +636 m s~^ (iLindegren fc Dravinsll2003l ). Therefore, the CB 
effect should be accounted for when determining a star's absolute RV, relative to its center 



of mass velocity (lAUende Prieto et al 



2OO2I : ILindegren fc Dravinsll2003l : lAUende Prieto et al. 



2009l : lde la Cruz Rodriguez et al.ll201ll ). and it should be taken into account when measuring 
the relative RV between t he two components of a binary star, when the required accuracy 
is better than 1000 m s"^ JPourbaix et al.ll2002[ l. 



Our goal here is to quantify the impact of the CB effect on Doppler shift measurement s 
during a planetary transit. This effect was already briefly mentioned by lWinn et al.l (120051 ). 
but its influence on the measured RVs was not studied. 



3. The Model 

We adopt a model where the planet is a spherical dark body, and the host star is a 
uniformly rotating spherical body with a limb darkened surface. We further assume that 
the host star has a uniform CB effect, where the radiation emanating from each position on 
its surface is blueshifted by Vcb directed perpendicular to the surface. A hidden assumption 
here is that the transiting planet size (or the stellar disk area the planet transverses during 
a single exposure) is much larger than the typical size of a granular cell, where for the Sun 
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the latter has a diameter of typically up to about 1,000 km (e.g., iRoudier fc Mullerl 119871 ) 



We note that the adopted model is somewhat simplified as it is meant to present a 
first step at understanding the impact of the CB effect on spectroscopic transit obser- 



vations. The model ignores pos sible more subtle effects (e.g., iLindegren fc Dravins 



2003 



de la Cruz Rodriguez et al.ll201ll ). such as meridional flows on t he stellar surfa ce (e.g.. lDravins 
19991 ). a possible different Vcp. for differe nt spectral lines (e.g.. lDravinslll982l ). or at different 
stellar l atitudes (iBeckers fc Tavlorlll980l ). or a dependence of the local observed blueshift on 
fi (e.g., Howard et al. 198o[ see also Section HJ. 



In our model the blueshift radial component of each surface element is fiVcB, so the 
absolute RV decrease induced by the CB effect is: 



Vo. 



CB 



Jo VbB/(/i)/W^C?/i 



CB 



where the integration is across the stellar hemisphere seen by the observer, /(/i) is the limb 
darkening law, and / is an order of unity coefficient. Since in our model Vcb is a constant 
we get: 



/ 



J^^I{n)fidfi 



(2) 



which equals 2/3 when there is no limb darkening. As the granulation pattern constantly 
changes, the value of / wil l vary with a time sc ale of a granulation cell lifetime, which for 
the Sun is about 10 min ( lAsplund et al.ll2000l ). However, even for short exposure times 
this variation is small if not negligible since it is averaged out by the integration across the 
stellar surface, containing many cells. Although granulation cell characteristics may change 
significantly for giant stars those are not targeted by high precision RV measurements. We 
note that the corresponding integral for the RM effect is zero when the planet is out-of- 
transit, meaning the origin of the RM effect — the host star rotation — does not affect the 
absolute RV, as expected. 

As an example, we have adopted an imaginary star-planet system whose parameters are 
typical for a transiting system with a close-in giant planet, and are listed in Tabled! We 
used a planetary radius and mass of Rp = 1.2 Rj and Mp = 0.8 Mj, and assumed a 3.5 
day orbit which is completely edge-on (i=90 deg), circular (e=0) and aligned (A=0 deg). 
For the host star w e assumed it is a GO type star with radius and mass of i?s = 1.10 Rq 
and Ms = 1.05 M© JCoxIboOoh . rotational velocity of Vrot = 2,000 m s"^ and CB velocity of 
Vcb = —500 m s~^. We calculated the Keplerian orbit of such a system, and at each phase 
during the transit we numerically integrated over the visible and eclipsed stellar surfaces in 
order to calculate both the RM and CB RV shifts. We assumed a quadratic limb darkening 
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law and determined the coefficients using the grids of I Clare tl (l2004f ) for the SDSS-g' band, 
where the host star emits most of its radiation. In this model the coefficient / in Eq. [1] is 
0.72 and Vq^cb = -360 m s~\ 

The resulting RV curves are plotted in Fig. [1] The top panel shows the RV curve during 
transit induced by the CB effect (blue), while the RM and Keplerian signals are ignored. 
The middle panel shows the CB RV curve (blue) and the RM RV curve (green) and their 
sum (black), while ignoring the orbital RV signal. For completeness we show in the bottom 
panel the CB and RM RV curves including the Keplerian orbit, with the same color code as 
the middle panel. 



4. Discussion 

During transit the planet blocks light coming from different regions of the host star's 
surface. This changes the integrated surface area in Eq. [T], and since surface elements with 
different sky-projected distance from the star center have a different RV component (following 
a different value for fi), the result of the integration varies, giving an anomalous RV signal 
during transit. This signal is superposed on the Keplerian and RM signals. We note that 
similarly to the RM effect, the RV variation originates from a variation in the spectral line 
centers following a variation in their shape, and not a shift of the lines as a whole, which 
is the case in orbital motion. Since the CB RV shift depends only on /i (and Vcb, which is 
constant) the CB RV curve is a symmetric function about the mid transit time (for circular 
orbits) . 

As shown by Fig. [T] top panel, the CB effect during transit induces a small anomalous 
RV signal, at the 1 m s~^ level, although it could be larger for other stars with larger Vcb- 
The signal reaches maximum at mid transit time since in our model the observed blueshift 
velocity is largest close to the center of the stellar disk, which is also where most of the light 
coming from the star is emitted due to limb darkening. 

The reason for the sign change in the CB RV curve can be understood as follows. During 
the transit start and end the planet blocks light coming from the stellar limb, where the radial 
blueshift velocity component is smaller than the average across the entire stellar disk. This 
in turn increases the overall blueshift, causing a decrease in the measured velocity. When 
the planet approaches the stellar disk center, the blueshift RV component of the echpsed 
stellar surface increases, thereby decreasing the average blueshift across the rest of the disk, 
inducing an increase in the average velocity in Eq. [H As seen in Fig. [1] top panel, the RV 
increase at mid transit is larger than the RV decrease close to transit start and end. This 
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reflects the geometry of tlie model and the limb darkened stellar surface. The maximum 
increase in RV can be approximated as: 



^CB « fVcB- ^ , (3) 

1 — 

where / is determined by the integrals in Eq. [2]and r is the planet to star radii ratio. For the 
example presented in Fig. [H Aqb = 2.3 m s~^. We caution that Eq. [3] is a rough estimate, 
and AcB decreases for inc lined orbits with i < 90 deg. An order of magnitude comparison 



to the RM RV amplitude (IGaudi fc WinnI 120071 . Eq. 6), gives: 



^CB VcB I . N 

-fi-RM Vrot 

This ratio is 25 % for the example presented in Fig. [H which uses Vrot = 2,000 m s~^, and 
it is larger for stars rotating more slowly. A more accurate estimate of the above ratio can 
be achieved by multiplying Vcb by / in the right hand side of Eq. HI resulting in a ratio of 
18 % for our adopted system, while the numerical integration gives a ratio of 12 %. 

The CB-induced RV curve does not depend on A or Vrot, meaning that on its own it 
does not hold information about the system's spin-orbit alignment. In fact, the CB RV 
signal should be present, in principle, even if the host star does not rotate at all. Since the 
measured RV curve includes both effects, the CB and RM, not accounting for the CB effect 
in the fitted model may result in a small bias of the spin-orbit angle. This is especially true 
for bright host stars that are also quiet, meaning show a low stellar activity level, as those 
stars allow a measured RV precision comparable to and even below the amplitude of the 
CB-induced RV signal. 

The symmetric shape of the RV curve induced by the CB effect distorts the shape of the 
RM signal. Specifically, for completely edge-on orbits (not necessarily aligned) or aligned 
orbits (not necessarily edge-on) the RM curve is an antisymmetric function about the mid 
transit time, so the combined RM + CB curve will have no definite symmetry, and the CB 
RV curve can not be completely absorbed into an RM model simply by changing Vrot or r. 
Therefore, the CB effect can bias the interpretation of the RV curve measured during transit 
and the estimate of A. For orbits close to edge-on the distortion is strongest during the 
central part of the transit where it increases the measured RV. Fig. [2] presents a zoomed-in 
view of the middle panel of Fig. [T] around the mid transit time. It shows that the CB + RM 
RV curve (black solid line) crosses the RV zero point (dashed black line), meaning it equals 
the Keplerian velocity, later than when the RM curve does, which is at the mid transit time 
(dotted vertical black line), when the RM effect briefly cancels out. This time difference is 
approximated by dividing the CB signal at mid transit time (Aqb) by the slope of the RM 
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RV curve: 

where r is the duration of the full transitl^, P is the orbital period, and Mg is the host star 
mass. Eq. O shows that the above time difference can easily be a few minutes or more. It is 
6 min for the system presented in Fig. [H Of course this time difference is dependent on the 
inclination angle, and for systems far from edge-on it will depend on A as well. However, 
for systems close to edge-on it will depend only weakly on A, unless A is close to 90 or 270 
deg. For example, for the system simulated here the above time difference is mimicked by 
the RM effect alone by increasing the impact parameter to about 0.25 and taking A ~ — 15 
deg. 

Therefore, an accurate estimate of the mid transit time and the expected Keplerian 
RV curve during transit are important for the correct interpretation of the spectroscopic 
transit RV curve and the identification of the CB effect. An accurate mid transit time can 
be obtained through simultaneous photometry — the measu rement of the photometric light 



curve of the same transit event observed spectroscopically (IHirano et al.l l2010l : IWinn et al. 



20 111 ) — since transit light curves are known to provide mid transit time estimates with an 
accuracy of 1 min and better. The required uncertainty on the Keplerian RV curve during 
transit is smaller than the CB effect amplitude. Such a low uncertainty on this RV curve can 
be achieved by obtaining RVs out-of-transit, both right before and right after the transit, 
and also on adjacent nights in orbital phases further away from the transit, to allow for an 
accurate orbital characterization, including specifically the RV slope during transit and the 
system's center of mass velocity. 

In addition, we note that for edge-on orbits with A of 90 or 270 deg the RM signal 
vanishes and the in-transit RV curve will be dominated by the CB effect. 

As we already mentioned the model we used here is probably a simplified version of 
the true behavior of the CB effect. Specifi cally, it was noted by several authors (e.g.. 



Beckers fc NelsonI 1 19781 : iHoward et al.l Il980l ) that the dependency of the observed local 



blueshift velocity scales as /i"VcB, where a > 1. This leads to a faster decrease in the 
observed Doppler shift towards the stellar limb, faster than in the model we adopted here 
where a = 1. The reason for this is the longer optical path in the stellar atmosphere closer 
to the limb. Using models with a > 1 in our calculations results in increasing Aq^. 



^The time from the second to third points of contact (jSeager fc Mallen-OrnelasI 
disk is completely enclosed by the stellar disk. 



20031 ) ■ when the planet's 
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Although beyond the scope of this paper, it is worth noting that the CB effect can 
influence the analysis of the RM effect in stellar binaries, where th e latter is measured 



through a careful analysis of the variation in the spectral lines shape (lAlbrecht et al.l 12011 
and references therein). 



5. Summary 

We have presented here the CB effect and its impact on the RV curve during a planetary 
transit. We adopted a simplified model where every surface element is radiating with the 
same blueshift velocity, directed along the normal to the surface. Using numerical integration 
across the observed stellar surface we calculated the RV curves for an imaginary transiting 
star-planet system, similar to some of the well known transiting systems with a close-in giant 
planet. 

The amplitude of the anomalous signal induced by the CB effect during transit is at 
the 1 m s~^ level for giant transiting planets. Not accounting for the CB effect can lead to 
inaccurate interpretation of the RV curve, possibly impacting the estimate of the spin-orbit 
angle, which will be especially visible in the case of bright and quiet host stars with accurate 
RVs. 

For close to edge-on orbits the CB RV signal during transit is strongest at mid transit 
time, where it causes a time shift between the mid transit time and the time when the 
RV curve crosses the Keplerian curve. Therefore, an accurate estimate of the mid transit 
time and a well determined orbit are important for an accurate analysis of the spectroscopic 
transit. 

So far the CB effect has been studied in detail only for the Sun. Future high-precision 
RVs will allow us to study it also for bright stars orbited by transiting planets, showing 
yet another way through which transiting planets can be used to probe their host stars. A 
measurement of the CB velocity for a large sample of stars will help to better understand 
processes such as convection and granulation, which in turn will help to better determine 
absolute RVs. 

We are grateful to Ake Nordlund for useful discussions, and to Dainis Dravins and 
Eric Agol for important comments. AS acknowledges support from NASA Grant Number 
NNX10AG02A. 
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Table 1. Table of parameters of the transiting planetary system used in the example 

presented in Figs. 1 and 2. 



Parameter 


Value 


Units 


Orbital parameters: 






Orbital Period, P 


3.5 


day 


Inclination, i 


90 


deg 


Eccentricity, e 







Spin-orbit angle, A 





deg 


Host star parameters: 






Rotational velocity, Vrot 


2,000 


m s"^ 


Convective blueshift, Vcb 


-500 


m s~^ 


Star radius, Rs 


1.10 


Rq 


Star mass, Ms 


1.05 


Me 


Planet parameters: 






Planet radius, Rp 


1.2 


Rj 


Planet mass, Mp 


0.8 


Mj 
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Fig. 1. — Top: The CB- induced RV curve during transit (blue), while the RM and Keplerian 
signals are ignored. Middle: The CB RV curve (blue) and RM RV curve (green) and their 
sum (black), while ignoring the Keplerian signal. Bottom: Similar to the middle panel, with 
the Keplerian orbit RV added to all three curves. All panels show RV in m s~^ versus orbital 
phase from mid transit time. 
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Time [min] 

Fig. 2. — Zoomed-in view of Fig. 1 middle panel around mid transit time. The different 
RV curves are represented using the same color code as in Fig. 1. The figure shows that 
for orbits close to edge-on or spin-orbit aligned, while the RM RV effect (solid green line) 
briefly cancels out and crosses the zero point (dashed black line, which is equivalent to the 
Keplerian signal) at mid transit time (vertical dotted black line), the RV curve including the 
RM and CB effects (solid black line) crosses the Keplerian signal at a later time. This time 
difference is quantified in Eq. [5] and equals 6 min in this example. The blue solid line is the 
RV curve induced by the CB effect. 



